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Abstract: 

Higgs-boson production in association with bottom quarks is an important discovery channel 
for supersymmetric Higgs particles at hadron colliders for large values of tan p. We present the 
complete 0{a) electroweak and 0{as) strong corrections to associated bottom-Higgs production 
through bb fusion in the MSSM and improve this next-to-leading-order prediction by known two- 
loop contributions to the Higgs self-energies, as provided by the program FeynHiggs. Choosing 
proper renormalization and input-parameter schemes, the bulk of the corrections (in particular 
the leading terms for large tan /3) can be absorbed into an improved Born approximation. The 
remaining non- universal corrections are typically of the order of a few per cent. Numerical results 
are discussed for the benchmark scenarios SPS lb and SPS 4. 
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1 Introduction 



The Higgs mechanism is a cornerstone of the Standard Model (SM) and its supersymmetric 
extensions. The masses of the fundamental particles, electroweak gauge bosons, leptons, and 
quarks, are generated by interactions with Higgs fields. The search for Higgs bosons is thus one 
of the most important endeavours in high-energy physics and is being pursued at the upgraded 
proton-antiproton collider Tevatron with a centre-of-mass (CM) energy of 1.96 TeV, followed in 
the near future by the proton-proton collider LHC with 14 TeV CM energy. 

Various channels can be exploited to search for Higgs bosons at hadron colliders. Higgs 
radiation off bottom quarks [1] 

pp/pp^bb0°+X (1.1) 

is the dominant Higgs-boson production mechanism in supersymmetric theories at large tan f3, 
where the bottom-Higgs Yukawa couplings are, in general, strongly enhanced. With cf)^ = H, h°, 
HO, AO we denote the SM Higgs boson or any of the neutral Higgs bosons of the Minimal Super- 
symmetric Standard Model (MSSM). Current searches for MSSM bottom-Higgs associated pro- 
duction at the Fermilab Tevatron widely exclude values tan /3 > 50 for light M^o « 100 GeV [2,3], 
depending in detail on the value of the Higgs mixing parameter /x. 

Two different formalisms have been employed to calculate the cross section for associated 
bb^o production. In a four-flavour number scheme (4FNS) with no b quarks in the initial 
state, the lowest-order QCD production processes are gluon-gluon fusion and quark-antiquark 
annihilation, gg — > bb 0° and qq — bb respectively. The inclusive cross section for gg — *■ bb 0° 
develops large logarithms ~ ln(/Lti?/mb), which arise from the splitting of gluons into nearly 
collinear bb pairs. The large scale /U^? ~ M^o corresponds to the upper limit of the collinear 
region up to which factorization is valid. The / m^j) terms can be summed to all orders 

in perturbation theory by introducing bottom parton densities. This defines the so-called five- 
flavour number scheme (5FNS) [4]. The use of bottom distribution functions is based on the 
approximation that the outgoing b quarks are at small transverse momentum. In this scheme, 
the leading-order (LO) process for the inclusive bb (p^ cross section is bb fusion, 

bb ^ 0° . (1.2) 

The next-to-leading order (NLO) cross section in the 5FNS includes 0{as) corrections to bb — *• 0° 
and the tree- level processes gb Ixp^ and gb hcfP. 

To all orders in perturbation theory the four- and five-flavour schemes are identical, but the 
way of ordering the perturbative expansion is different, and the results do not match exactly at 
finite order. However, numerical comparisons between calculations of inclusive Higgs production 
in the two schemes [5-8] show that the two approaches agree within their respective uncertainties, 
once higher-order QCD corrections are taken into account. 

There has been considerable progress recently in improving the cross-section predictions for 
inclusive associated bb^^ production by calculating NLO-QCD [5,7] and SUSY-QCD [9] correc- 
tions in the four- flavour scheme, and NNLO QCD corrections [10,11] in the five- flavour scheme. 
The inclusion of higher-order effects is crucial for an accurate theoretical prediction and, eventu- 
ally, a determination of Higgs-boson parameters from the comparison of theory and experiment. 
In this paper we further improve the cross-section prediction and present the first calculation 
of the complete 0{a) electroweak corrections to associated bottom-Higgs production through 
bb -> / in the MSSM. 

The complete one- loop QCD and electroweak corrections for the decay of MSSM Higgs bosons 
to bottom quarks have been presented in Ref. [12] more than a decade ago. While the predictions 
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for Higgs-boson decays have been improved and refined in recent years, supersymmetric QCD 
and electroweak corrections to tlie production cross section have so far only been investigated at 
the level of universal corrections for large tan/3 (see e.g. Refs. [3,13]). 

In Section [2. II we shall set the notation for the supersymmetric model. The calculation of the 
NLO QCD and electroweak corrections is described in some detail in Sections 12.21 12.31 and 12. 4[ 
Numerical results for MSSM Higgs-boson production at the LHC are presented in Section [3l We 
conclude in Section [H The Appendices provide details on the scenarios of the supersymmetric 
model under consideration and present further numerical results. 



2 Radiative corrections 



2.1 Tree-level Yukawa couplings and cross section 

In the Standard Model, Higgs production from bottom-quark fusion is governed by the interaction 
term 

>CbbH = -A^''bbH, (2.1) 

where Ab is the bottom-quark Yukawa coupling and H denotes the field of the physical Higgs 
boson. The corresponding mass term for the bottom quark is generated by the Higgs vacuum 
expectation value v, leading to the tree- level relation 

where e is the electromagnetic coupling constant, Sw the sine of the weak mixing angle, and Mw 
the W-boson mass. 

For the MSSM we will follow the conventions of Ref. [14], where the two Higgs doublets are 
denoted as 

\^^{v,+h«+ix'i)) \ -h„- ) 

The bottom quarks couple to Hd, giving mass to the down- type quarks via the vacuum expec- 
tation value Vd- Masses for up- type quarks are generated by a second Higgs doublet Hu with 
vacuum expectation value Vu_- Considering the MSSM without CP-violating phases, the CP-even 
neutral Higgs-boson fields h*^ and H'^ are linear combinations of h[} and h[}. One conventionally 
defines the Higgs mixing angle a by writing 

/" — c \ / ViO \ 

(2.4) 

Here and in the following, we will frequently use the notation Sa = sin a, Cq, = cos a and 
generalizations thereof. The vacuum expectation values are parameterized according to = 
V sin/3, Vd = V cos/3, i.e. 

ti3 = tan /3 = — and v'^ = v'^ + v'^ . (2-5) 

Taking the pseudoscalar Higgs mass M^^o and t/j as input parameters of the MSSM Higgs sector, 
one finds the tree-level relation 

_ M2 + M| 
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with S2a < 0. The fields of the physical neutral pseudoscalar Higgs boson A*^ and the neutral 
would-be Goldstone boson G*^ are given by 





(2.7) 



C/3 -SfS 
8/3 C/3 

Consequently, the bbh*^, bbH*^, and bbA*^ couplings read 

\hO _ -Sg mb _ _\SM Sq 

\H0 _ Cg TTT-b _ ^SM Cg (2 a^ 

'^b - -^U^ - \ c^^ 
pP _ -Sf3 mb _ ,sM^ 

Hence, the Yukawa couplings are enhanced for large values of tp. Note that for large masses of 
the pseudoscalar Higgs boson, h*^ is known to be SM like and Sg — > — C/3. 
The leading-order partonic cross sections are given by 



c2/c2 for hO production, (2.9) 
t} AO 



where 4>° = (h°,H°,A°), r = M|o/s, is the partonic CM energy, and the incoming bottom 
quarks are treated as massless particles in accordance with QCD factorization. 

2.2 SUSY-QCD corrections 

The QCD corrections to the process bb h*',H'^,A'^ are known to NNLO [10,11], with a small 
residual QCD factorization and renormalization scale uncertainty of less than 10%. If one 
chooses the renormalization and factorization scales as = Mh and = Myi/i, respectively, 
the impact of the genuine NNLO QCD corrections is typically less than 5% [11]. We have 
reproduced the NLO QCD result and extend previous analyses by including the 0{as) SUSY- 
QCD corrections from virtual squark and gluino exchange. 

The MS scheme has been adopted for the renormalization of the bottom-quark mass mb and 
for the factorization of initial-state collinear singularities. The renormalization of the bottom- 
Higgs Yukawa coupling is fixed in terms of the bottom-mass renormalization. In order to sum 
large logarithmic corrections oc ln(mb/^_R) we evaluate the Yukawa coupling with the running 
b-quark mass mb(^_R) [15]. 

The 0(0^) SUSY-QCD corrections comprise self-energy and vertex diagrams induced by 
virtual sbottom and gluino exchange, as shown in Fig. [H It is well known [16-19] that the 
SUSY-QCD corrections are enhanced for large t^. This important effect can be qualitatively 
understood as follows. Unlike down-type quarks, which only couple to the down-type Higgs 
field at tree level, the down-type squarks also couple to the up-type Higgs field via terms in the 
superpotential. The corresponding coupling strength is proportional to the enhanced Yukawa 
coupling mb/vd times the Higgs mixing parameter fi. Hence, the vacuum expectation value 
of Hu leads to a mixing term in the sbottom mass matrix proportional to /it/j, which dominates 
the sbottom mixing angle Gj^, 

sm(2eg) = 2 2 ' (2-10) 

m- — m- 

bi b2 
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(a) 



(b) 



Figure 1: Self-energy (a) and vertex (b) corrections from gluino exchange. Sbottom mass 
eigenstates are denoted as with i = 1,2. 



where denotes the soft-supersymmetry-breaking trihnear scalar coupling and bi^2 are the 
sbottom mass eigenstates. The factor sin(2Gg) directly enters the (scalar part of the) b-quark 
self-energy, which in turn enters the Yukawa coupling renormalization via the b-mass counterterm 
5mb. The corresponding mass shift is usually denoted by — mbAb with 

^h = ^^rn~gnti3l{m^^,m^^,m~g), (2.11) 

Cp = 4/3, and the auxiliary function 

Here, rrig is the gluino mass. As shown in Ref. [20] by power counting in agt^, the contribution 
of Ab can be summed by the replacement 

^ (2.13) 

1 + Ab ^ ' 



in the bottom Yukawa coupling. As explained above, the loop-induced coupling of the up-type 
Higgs to bottom quarks also involves a factor astp. The full contribution to the ^'^bb vertex 
receives an additional factor {cq, s^, c^} from the Hu part in (fp = {h°, H°, A°}. Power counting 
in tp shows [20] that the t^-enhanccd vertex corrections of the form (a^. t^)" are one-loop 
exact, i.e. they do not appear at higher orders (n > 2). Collecting all t/3-enhanced corrections to 
the bottom Yukawa couplings leads to the effective couplings [20,21], where 



'>% _ _Sa 1- Ab/(t/3ia) 

(2.14) 



ASM 1 + Ab 

Ag° ^ Cg 1 + Ab ta,/tp 

ASM 1 + Ab 

1 - Ab/t| 
ASM - i + Ab • 

Note that A{J° is still SM like for large M^o, independent of the large-t^ summation owing to 
t/sta ^ — 1 in this limit. The summation formalism can be extended [21] to include corrections 
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proportional to the trilinear coupling in (|2.10p . However, these corrections turn out to be smah 
and summation effects may thus safely be neglected for the MSSM scenarios under consideration 
in this work. 

To combine the features of the above effective treatment with the complete one-loop SUSY- 
QCD calculation, we modify the renormalization scheme to absorb the above corrections into a 
redefinition of the bottom mass in the Yukawa coupling. Hence, an additional counterterm 

is added for h^, H", and A'^ production, respectively, to remove the t/3-enhanced contributions 
from the explicit one-loop result in order to avoid double-counting. We use the convention of 
Ref. [12] where mb,o = "T-b + Sm\^. 

As we shall demonstrate in the numerical analysis presented in Section [3l the SUSY-QCD 
radiative corrections are indeed sizeable at large t^. After summation of the t^-enhanced terms, 
however, the remaining one-loop SUSY-QCD corrections are negligibly small, at the level of per 
mille and below. 



2.3 Electroweak SM corrections and calculational details 

The electroweak corrections naturally decompose into a purely photonic, QED-like part 5qed 
and the remaining weak contributions ^weak^ ^ew = Sqed + <^weak- Each of these contributions 
forms a gauge-invariant subset of the C'(a)-corrected cross section. The photonic corrections 
due to virtual photon exchange and real photon emission can be obtained from the QCD results 
by appropriately adjusting colour factors and electric charges. For the QED renormalization of 
the bottom mass we use the on-shell scheme, because electroweak running effects beyond the 
one-loop level are negligible. 

The divergences due to collinear photon emission from the massless b-quarks are removed by 
mass factorization as in QCD, i.e. by a redefinition of the bottom parton densities according to 

where A = 1/e — 'Je + ln(47r) is the standard divergence in L> = 4 — 2e dimensions, "Je is Euler's 
constant, Qb = ~l/3 is the electric b charge, and he denotes the QED factorization scale which is 
identified with the QCD factorization scale but is chosen independently of the scale /x introduced 
by dimensional regularization. The factor 3 in the second line stems from the splitting of the 
photon into bb-pairs of different colour. Furthermore, 

P<i<i(.^) = \^ and Pg^iz) = + {I - zf (2.17) 



-C„iz) 

Cq^{z) 



(2.16) 
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are the quark and photon sphtting functions, respectively, and Cqq, Cqj the coefficient functions 
specifying the factorization scheme. Following standard QCD terminology one distinguishes MS 
and DIS schemes defined by 



„ , , , 1 - z 3\ 9 + 5z 
P,,(z) In^-- + 



z 4/ 4 



(2.18) 



1 - z 



C^^^{z) = Pqy{z) In 8z^ + 8z-l . 



The equivalent factorization procedure using fermion masses as regulators instead of dimen- 
sional regularization can, e.g., be found in Ref. [22]. In our numerical analysis we employ the 
MRST2004qed parton distribution functions [23], which include 0{a) corrections defined in the 
DIS factorization scheme [22]. The MRST2004qed parameterization also provides a photon den- 
sity necessary to compute the hadronic cross section for the 0{a) photon- induced processes 
'yh ^ (p^h and ^ (p^ h. 

We have calculated the NLO QCD and QED corrections using dimensional regularization 
and alternatively using a mass regularization for the collinear divergences, where we employed 
the methods of Refs. [22, 24] for the mass regulators. 

Similar to the QCD case, the QED corrections are universal for the production of SM and 
MSSM Higgs bosons. Their size is quite small since the potentially large correction due to 
collinear photon emission is absorbed into the parton distribution function f\j{x,fj,F). We shall 
discuss numerical results in Section [3l 

Let us now turn to the remaining electroweak corrections 5weak- As we shall discuss in more 
detail in Section 13.41 the small but finite bottom mass can induce sizeable corrections in the 
electroweak sector of the MSSM due to additional (t/j-enhanced) bottom Yukawa couplings in 
loops. Hence, while we neglect the b-mass at tree level and for the QCD and QED corrections 
as required by QCD factorization, we keep the finite bottom mass rrib in the calculation of the 
relative one-loop weak correction (5weak- Thus, our result for J^eak also contains kinematical rrib 
effects that formally lie beyond the accuracy of the 5FNS-calculation. However, these effects are 
small. Different choices for the numerical value of the bottom mass used within the calculation 
of the relative one-loop correction lead to results which formally differ by NNLO effects. Because 
the rrib-dependence is dominated by the strength of the Yukawa coupling, we have chosen the 
running bottom mass (as defined after summation of t/j-enhanced terms) as input. 

Both in the SM and in the MSSM, the one-loop electroweak corrections have been calculated 
independently using two different approaches: In one approach, the Feynman-diagrammatic 
expressions for all self-energies and one-loop vertex diagrams have been generated using the 
program package FeynArts [25]. The calculations have then been performed with the help of 
the program package FormCalc [26], and the loop integrals have been evaluated numerically with 
LoopTools [26]. In a second approach, also starting from the amplitudes generated by FeynArts, 
all calculations, including the evaluation of the loop integrals, have been performed with a com- 
pletely independent set of in-house routines. The two calculations are in mutual agreement. We 
note that the regularization of the complete electroweak MSSM corrections has been performed 
using both constrained differential renormalization as implemented in FormCalc as well as di- 
mensional reduction. Both regularization procedures are known to be equivalent at the one-loop 
level [26], and the results of the two calculations are in agreement. We refrain from displaying the 
complete analytic results and restrict ourselves to a discussion of the renormalization conditions 
and the input-parameter schemes. 
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Using standard notation, the vertex counterterm at one-loop order is given by 
.bbH _ , 5Zh 6m^ bZ\ + bZ\ 6s^ 6M^ 



Employing the on-shell renormalization scheme, this results in 



(2.20) 



where denotes the scalar and vector part of the b-quark self-energy, respectively, and Tl 

refers to the derivative of the self-energy with respect to the external momentum squared. The 
relation = 1 — M-^/M| is used to determine (5s w Here and in Section [2^ we only consider 
the real part of the self-energies and follow the conventions of Ref. [27]. For compactness, the 
transverse parts of gauge-boson self-energies are simply written as Ew, Ez, etc. 

The different input-parameter schemes are specified by the choice of bZ^.. In the a(0)-scheme, 
using the low-energy fine-structure constant a(0) as input, we have [280 



*Z,U,„,.^E;,0) + ^^. (2.21, 



In the a(0)-scheme, bZf, contains logarithms of the light fermion masses inducing large corrections 
of the form aln(mj/s), which are related to the running of the electromagnetic coupling a{ff) 
from Q = to a high energy scale. In order to correctly reproduce the non-perturbative hadronic 
part in this running, which enters via n^(M|) - n^(0) with n^(A;2) = S-y(/c2)/A:2 denoting the 
photonic vacuum polarization, we adjust the quark masses to the asymptotic tail of n^(A;2). In 
the a(Mz)-scheme, using a(Mz) as defined in Ref. [29] as input, this adjustment is implicitly 
incorporated, and the counterterm reads 

^^e|,(Mz) = '^^e|„(o) " ^Aa(M|) , (2.22) 

where 

Aa(Q2) = n/^*(0) - Ren/^*(Q2) , (2.23) 

and 11;^^ denotes the photonic vacuum polarization induced by all fermions other than the top 
quark. Hence, in the a(Mz)-scheme, the final result does not depend on the above logarithms of 
the light fermion masses. In the G^-scheme, a is determined from the muon-decay constant 
according to 

= ^G,Mls\ ^ ^(0)(1 + ^,)_ (2.24) 

TT 

The radiative QED corrections for muon decay in the framework of the effective 4-fermion inter- 
action are already encoded in the numerical value for G^. The additional corrections from a full 
one-loop SM calculation are taken into account through Ar [30] according to 

5Ze|G^ = 5Ze|,(0)-^Ar. (2.25) 



^We follow the convention of Haber and Kane [14] to define the covariant derivative for SU(2)l, i.e. the sign in 
front of E^z in (p:2T1l differs from Ref. [28] . 
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Since Aa(M|) is explicitly contained in Ar, the large fermion-mass logarithms are also absent 
in the G^-scheme. Moreover, since the lowest-order cross section is proportional to a/ for the 
production of all the three Higgs bosons, in G^-parameterization it absorbs the large universal 
correction Ap from the /O-parameter, which is oc G^m^ and represents a part of Ar. Dividing 
Eq. (j2.24p by s^, one easily sees that ac^/syi^ absorbs Ar and thus also Ap. We will use the 
G^-scheme unless stated otherwise. We have also performed two independent calculations for 
Ar in the MSSM, using either constrained differential renormalization or dimensional reduction, 
and find agreement with the result of Ref. [31] . 

In the SM, the Higgs mass can be defined by an on-shell renormalization condition. The 
wave-function renormalization of the Higgs-boson field is conveniently chosen in the on-shell 
scheme, 

5Zh = -Sij(M2), (2.26) 
where Sh is the Higgs-boson self-energy. 



2.4 Electroweak MSSM corrections 

The photonic corrections in the MSSM do not change with respect to the SM case. For the 
conventions and for the renormalization of the MSSM Higgs sector, which is more involved, we 
essentially follow Ref. [27] H In particular, a proper renormalization scheme has to be specified 
to determine the vertex counterterm S^^, cf. (j2.19p . including the renormalization of tp. The 
wave-function renormalization for the Higgs doublet fields is usually defined by 

H, ^ Z^f H, = H, (^1 + i 6Zn,^ (2.27) 

with i = u, d. For the vacuum expectation values of the Higgs fields one defines 

Vi ^ Zl/^{vi - 5vi) = Vi(l + \ 5Z^^ - ^) , (2.28) 

where the last equation is valid to one-loop order with Zh . = 1 + ^Z-^. . The freedom of wave- 
function renormalization can then be used to impose the condition 

5vvL 6va 

— = — (2.29 

leading to 

^ = i (5Zh. - 5Zh J . (2.30) 
Hence, for the MSSM Higgs sector the counterterm 5Zy{ in (j2.19p reads 

5Zn = i 5Zh, + 4 ^ = ^ (4 5Zn, + 4 <5Zh J . (2.31) 

Note, that (5Zh includes all parts of the vertex counterterm that are related to the Higgs sector, 
i.e. 5tp as well as the wave-function renormalization counterterm SZn^. The quantity (5Zh should 



^For clarity, we specify our conventions for some field-theoretic quantities where the conventions of Ref. [27] 
might be unclear. Explicit tadpole vertex functions for Higgs fields are denoted as T** = iT''' , i.e. T"^ differs 
from Ref. [27] by a global sign. The A°Z mixing self-energy is derived from the vertex function ^(fc, — fc) — 
fcfiEA_oz(fc^), where k is the incoming A" momentum. 
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not be confused with the wave-function renormahzation constants for the physical Higgs fields. 
Since the counterterm dZu is universal for the production of h", H'^, and A'' via down- type 
quarks, the whole vertex counterterm ()2.19p is also universal. 

We consider two renormahzation schemes for t^: 

i) Following Dabelstein [27] and Chankowski et al. [32], a vanishing on-shell A^Z-mixing can 
be used renormahzation condition, i.e. 

Saoz(M1o)=0, (2.32) 

where 

Eaoz(A:') = $]aoz(A;2)-MzS2/3^ (2.33) 

1/3 

denotes the real part of a renormalized self-energy defined according to the conventions 
in Ref. [27]. To fix the second wave-function renormahzation constant, one demands the 
on-shell condition 

±'j^o{Mlo) = (2.34) 
for the residue of the A^-boson propagator. From the last two equations, one finds 

6Zn^ = -T.'^.{Ml,) + ^ Saoz(mIo) , 

1 



(2.35) 



5Zn, = -S'Ao(Mio) - ^ Saoz(MIo) , 

which defines the DCPR scheme [27,32] for the tp renormahzation [see (|2.30p ]. 

Alternatively, in the DR scheme, the counterterm for is proportional toA = l/e — 7£;-|- 
ln(47r). Hence, one can convert ()2.35p to the DR scheme by setting the remaining finite 
pieces of Sa^'z to zero. Accordingly, is a renormalization-scale-dependent quantity, i.e. 
the input for a given model has to be fixed at a given scale /i/j. To one-loop order, the 
conversion of the tp input parameters from the DR scheme to t^^^^ in the DCPR scheme 
is given by 

DCPR 



+ o;i^_„2 ^AOz(^^Ao] 



(2.36) 



1 

2Mzc;3 

where S^"^ denotes the finite pieces of the self-energy in the DR scheme. In this work, we 
always use the renormahzation condition (j2.34p . also if we use DR to renormalize tp. 

In analogy, there are the DCPR and the DR schemes for the renormahzation of the mass 
Mpfi of the pseudoscalar Higgs boson. The DCPR scheme uses the on-shell renormahzation 
condition 

Sao(A^1o) = 0, (2.37) 

while the DR scheme is again defined by setting S^o(-^ao) *° ^^'^^ ™ mass counterterm for 
the pseudoscalar Higgs boson. In this work, we use the on-shell scheme for the renormahzation 
of A'^. Supersymmetric models are usually defined in terms of DR parameters [33]. Hence, we 
have to calculate the on-shell A*^ mass from the corresponding scale-dependent DR parameter. 
For a given parameter set, we determine M^o from the zero of the inverse A" propagatoiH 

{M°^,f - {M^f + T.%{{M°^of ) = (2.38) 



^In contrast to Ref. [19] , but in line with Ref. [33] , we assume that all tadpole contributions to the mass of A" 

fDR 
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which corresponds to a given DR mass for AO. Here, T.^o{k'^) is calculated from an MSSM 
parameter set using M^q as input. We start with M^o = for the self-energy calculation and 
iterate until self-consistency is reached. If is renormalized in the DCPR scheme, (|2.36|) is used 
to also find t^^^^ self-consistently along with M^q • The one-loop shift of the numerical value of 
M^o is particularly important because it enters already at tree level through the determination 
of the mixing angle [see (|2.6p ]. 

For completeness, we also state the remaining renormalization conditions for the tadpoles. 



TjjO -|- (5tho = and Tjjo -|- (5iH0 



0. 



(2.39) 



which ensure that t>u and correctly minimize the one-loop potential. As in the SM, the masses 
for the W and Z boson are renormalized by on-shell conditions. 

Including the corrections to the Higgs external legs, the partonic cross section to one- loop 
order is given by 



(1) 



a 



(1) 
HO 

(1) 
AO 



_tree 



„tree 
0-HO 



Zho 



1 - ^hORO - 



^HOhO — 

Cn 



+ 2Re[ 1 -ZhOHO- 



A 



hO 



+ 2Re{ I-Znoho-^ ) Aho 



(2.40) 



a^T-ZAo(l + 2 ReA 



AO; 



where A^o denotes the relative one-loop vertex corrections including the corresponding coun- 
terterms. Depending on the renormalization scheme, A^o also includes the contributions from 
ZAO -mixing and G'^A'^-mixing. The Z factors are given by [12] 



V2 



{fc2) 



Zho 



Zao 



1 



^HOhO 



/c2-m2o+Ei^o{A:2) 



(2.41) 



1. 



i + s^«(fc2) 

The mixing of the CP-even Higgs bosons is determined by 

ZhOHO = 



^HOhO 



SHOh0(M-^O; 



where m^o denotes the tree-level masses and M^o the one-loop masses 
zeros of the inverse one-loop propagator matrix determinant 



(2.42) 



h°, HO), i.e. the 



+ Sho(A;") A:^ -m^o + SHo(r) -S 



^HOhO \ 



0. 



(2.43) 
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The renormalized Higgs self-energies in turn are given by 

Sho(^') = SHo(fc') + {5ZuA + SZuA) - S^lo , (2.44) 
5^HOho(^^) = ^uoh»{k'^) + k"^ SaCa {SZ^^ - dZu^) - 5mHOho , 
where the Higgs mass counterterms reac^ 

r 2 2 ^ fj,j2 \ , 2 X- n\j2\ ^^l-a^P-a ^ S/3_q, (1 + c|,_^) 

6m^o = c^_« Sao(Mao) + s^+„Sz(Mz) - Tho + 

- (M^o c2_^ + M| S:,o(Mlo) + Mz ^^^^^^ Saoz(MIo) , (2.45) 

= s,.^ Sao(Mao) + c,+„Sz(Mz) + T^HO - 2Mz.wCw 

- (M^o + M| S:,o(M1o) - Mz £?f^±^ Saoz(MIo) , (2.46) 

2p 

3 3 

+^ (Mio S2(p.a) + Mi S2(p+a)) S'ao(M1„) + Mz ^ Saoz(MIo) . (2.47) 

If is renormahzed in the DR scheme the finite parts of S^oz have to be set to zero in the above 
formulas. Note also that M^o is the on-shell A^-mass in these formulas, because we translate 
the M^o input always to the on-shell scheme before the actual loop calculation. As mentioned 
before, in all the above equations we only consider the real parts of the self-energies. Neglecting 
the imaginary part does not spoil the one-loop accuracy of our calculation, since the imaginary 
parts formally only enter at higher orders H 

In complete analogy to SUSY-QCD in Section 12.21 there are t/j-enhanced corrections in the 
electroweak sector which can be numerically sizeable and which should be summed [20] to all 
orders. Higgsino-stop loops lead to a contribution to Ab of the form 

2 

= o"7^2 /i I{m^, , m-,^ , /i) , (2.48) 
where m? denote the masses of the stop mass eigenstates, and is the usual trilinear soft 



breaking parameter. From wino-higgsino-stop and wino-higgsino-sbottom loops we find 



= - ^^^2 /it/3 2 cos^ e,- I{mi^ , /i, M2) + 2 sin^ 6,- /(m,-^ , /x, M2) 



+ cos' I{m^^ , /X, M2) + sin' /(m^^ , /i, M2) 



Up to some (known) sign errors in Srn^a^o, the Higgs mass counterterms can be also found in Ref. [27]. 

^While an inclusion of such imaginary parts in mass determinations of unstable particles is straightforward, 
a consistent inclusion of such width effects in the description of particle reactions requires an inspection of full 
resonance processes including production and decay of the unstable particles. This is beyond the aimed level of 
precision of this work. 
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where M2 is the soft SUSY-breaking mass parameter for the wino, denotes the stop mixing 
angle, and the auxihary function I{a,b,c) is given in (I2.12|) . Finally, bino loops contribute 



A? 



■j^Mi tp [3(1 + sin^ G^) I{m^^ , Mi) + 3(1 + cos" 6^) I{m^^ , /x, Mi) 
+ 2/(m^^,m^^,Mi)j , 



(2.50) 



where Mi is the soft SUSY-breaking mass parameter for the bino. The full result for Ab is then 
given by 

Ab = A| + Ar"'' = A| + Ap + Af + A^, (2.51) 



where A^ denotes the gluino contribution of (j2.1ip . We have recalculated these results and 
find agreement with Ref. [20]. In contrast to Ref. [20] we also include the bino terms in the 
summation. However, they are indeed numerically small. 

To further improve the accuracy of the calculation, we also include two-loop contributions 
to the self-energies as provided by the program package FeynHiggs [35] (version 2.3.2). As for 
the one-loop part of the self-energies we find perfect agreement between our calculation and the 
results obtained with FeynHiggs 

In addition to using a properly defined b-mass, one often absorbs parts of the radiative 
correction related to the Higgs external leg using an effective mixing angle aes = a + 6a at 
tree level. Here, we follow FeynHiggs and define 5a to be the angle which diagonalizes the 
loop-corrected Higgs mass matrix, i.e. 



Sho (m^ 



+ <o)/2) 



m 



HO 



ho + m^o)/2) 
SH"("iHo) 





( Ml, 






\ 





(2.52) 



where the renormalized Higgs self-energies are evaluated at the given momenta. We define an 
improved Born approximation ctiba which includes the leading higher-order corrections through 
the running b-mass, the summation of t/3-enhanced terms and the effective mixing angle. 



C^IBA 



V27rG^mb(/ij;) ^ 
6M2 



6{1 




1 + Ab 
1 + Abt^,Jt^ 

AbA|^ 



hO 

for H" production, 
AO 



(2.53) 



where (Tiba denotes the partonic cross section, cjp = (hO,HO,AO), r = M'^^/s^ and \/| is the 
partonic CMS energy. Here mb(/^_R) is the running MS bottom mass in QCD, Ab comprises 
the t^-enhanced terms of the supersymmetric QCD and weak corrections as specified in (j2.5ip . 
and Oeff is the effective mixing angle defined in (j2.52p . We will compare our full result to this 
approximation in Section [3l 



''Since FeynHiggs does not directly support our renormalization scheme ii) (the Higgs field renormalization is 
different), we had to implement this scheme in FeynHiggs ourselves. For the leading two-loop corrections included 
in FeynHiggs, neither tan/3 nor the Higgs fields are renormalized at the two-loop level. Thus, it is consistent to 
add the FeynHiggs two-loop correction to the one-loop result in our renormalization schemes. 
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3 Phenomenological analysis 



3.1 SM input parameters 

For our numerical predictions, we essentially use the SM input parameters [36] 



a 


= 1/137.03599911, 


q(Mz) 


= 1/128.952, 


G, 


= 1.16637 X 10"^ 


Mw 


= 80.403 GeV, 


Mz 


= 91.1876 GeV, 






me 


= 0.51099892 MeV, 




= 105.658369 MeV, 


mr 


= 1.77699 GeV, 


mu 


= 66 MeV, 


rric 


= 1.2 GeV, 


mt 


= 174.2 GeV, 


md 


= 66 MeV, 




= 0.15 GeV, 


mb(rnb) 


= 4.2 GeV. 



_ (3.1) 

Here, mb(mb) is the QCD-MS mass for the bottom quark while the top mass nit is understood as 
an on-shell mass. For the QED renormalization of the bottom mass we use the on-shell scheme, 
as mentioned above, with an on-shell mass mb = 4.53 GeV calculated from mb(mb) to one- loop 
order in QCD. The masses of the light quarks are adjusted such as to reproduce the hadronic 
contribution to the photonic vacuum polarization leading to a{Mz) of Ref. [37] . They are relevant 
only for the evaluation of the charge renormalization constant 6Ze in the a(0)-scheme. The CKM 
matrix has been set to the unit matrix. For the calculation of the hadronic cross sections we 
have adopted the MRST2004qed parton distribution functions [23] at NLO QCD and NLO QED, 
with the corresponding as{Mz) = 0.11899. The top quark, squarks, and gluinos are decoupled 
from the running of the strong coupling asinR). We choose the renormalization and factorization 
scales as /i/j = Mh and fip = Mn/A, respectively. As mentioned above, with this specific choice 
the QCD NNLO radiative corrections are at the percent level [11]. 

3.2 The SM cross section 

The NLO cross section predictions for associated bbH production at the LHC, including QCD, 
QED and weak corrections, are shown in Table [U for the three different input-parameter schemes 
introduced in Section [2.31 The QED corrections are very small after potentially large contribu- 
tions from collinear photon emission have been removed by mass factorization as described in 
Section 12.31 As expected, the inclusion of the electroweak corrections reduces the scheme de- 
pendence. In the following we will adopt the G^-scheme, where the value of the electromagnetic 
coupling is derived from muon decay according to ac,, and where the influence of the light-quark 
masses is negligible. Note that the relative 0{a) QED corrections are evaluated with a(0), irre- 
spective of the chosen input-parameter scheme, because the relevant scale for the bremsstrahlung 
process is set by the vanishing virtuality of the emitted real photon. 

3.3 MSSM input 

For our numerical analysis, we will focus on the benchmark scenario SPS 4 [38] which is charac- 
terized by a large value of tp = 50 and a correspondingly large associated production cross section 
bb — > H^jA" at the LHC. Results for the alternative benchmark scenario SPS lb with tp = 30 
will be presented in Appendix [Bj Both the SPS 4 and SPS lb input parameters are specified 
in Appendix El Note that searches at the Fermilab Tevatron are not sensitive to scenarios with 
H*^ , A'^-masses as large as the masses in SPS 4 or SPS lb. 

From the SPS DR input parameters we calculate M^o a-iid, if we work in the DCPR scheme, 
also ifl"^^^, as specified in Section [2.41 The corresponding renormalization scale ^r(DR) used 
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a[ph] 


LO 


QCD 


QCD+QED 


full SM 


a(0)-scheme 
G^-scheme 
a(Mz)-scheme 


0.02309 
0.02390 
0.02453 


0.02868 
0.02968 
0.03047 


0.02863 
0.02963 
0.03042 


0.02901 (25.6%) 
0.02924 (22.3%) 
0.02932 (19.5%) 



Table 1: The LO and NLO SM cross section pp — > (bb)H + X for a Higgs boson with 
Mh = 300 GeV at the LHC {^/s = 14 TeV). Results are presented for the three different 
input-parameter schemes defined in Section 12.31 The MRST2004qed parton distribution 
function and NNLO-QCD running for the b-mass have been adopted for the NLO as well 
as the LO cross sections, and the renormalization and factorization scales have been set 
to /XR = Mh and //_f = Mh/4, respectively. "QCD" denotes the NLO QCD corrections 
only, "QCD+QED" also includes photon exchange and emission as well as the initial state 
containing a photon. The "full SM" prediction includes all 0{as) and 0{a) corrections; 
the corresponding relative correction is indicated in parentheses. 



in the electroweak part of the calculation is specified by the SPS scenario (see Appendix . 
Note, that it differs in general from the renormalization scale in the QCD part of the calculation 
which is set to the mass of the produced Higgs boson. The MSSM tree-level relations are used to 
determine the sfermion and gaugino masses and mixing angles that enter the one-loop corrections. 
The Higgs masses and the Higgs sector mixing angle a [see ()2.6p ] which enter the calculation of 
the loop diagrams are also obtained according to the tree-level relations. 

While the bottom mass that enters the Yukawa coupling at LO is fixed by the requirement to 
account for dominant NLO corrections, different definitions for the bottom mass in the relative 
NLO correction change the result only beyond NLO. As argued in Section 12.31 we shall use 
the running bottom mass as defined after summation of t^-enhanced terms as inpu10. The 
running mass is needed as an input for the determination of the MSSM parameters, but it also 
depends on these parameters through the QCD renormalization scale fiR = M^o and through 
the t^-enhanced corrections. The b-mass is thus calculated using an iterative procedure, starting 
from some initial guess for its value, and using the resulting b-mass as input for a refined 
determination of the MSSM parameters, until self-consistency is reached. Note, however, that 
this running mass depends on the process under consideration through the choice of scale and 
through the i^-enhanced corrections. In order to avoid the proliferation of input masses, we adopt 
the running b-mass associated with H*^ production in the relative corrections to all processes, 
e.g. TTib = 2.24 GeV for SPS 4 in the DR scheme. 

We use the two-loop-improved Higgs masses for the kinematics, e.g. in the tree-level cross 
section or when they appear as external momenta in the on-shell vertex-correction. Using the two- 
loop self-energies from FeynHiggs [35], these two-loop improved on-shell Higgs-boson masses in 
the SPS 4 scenario are given by M^o = 115.66 GeV, Mho = 397.72 GeV, and Mao = 397.67 GeV 
for the renormalization scheme ii) introduced in Section 12.41 i.e. what we call DR scheme. (As 
described in Section 12.31 we only consider the real part of the self-energies when we determine 



^While the electroweak corrections are calculated using dimensional reduction, we do not convert the running 
MS bottom mass to the corresponding DR mass to define the input value for the relative one-loop correction. The 
difference is of higher order. 
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hO 


H° 


A° 


0"IBA[pb] 


0.65 


15.39 


15.40 


^qcd[%] 


36.25 


21.48 


21.48 


Sqed[%] 


-0.13 


-0.23 


-0.23 


<5mssm-qcd[%] 


-0.03 


0.08 


0.07 


5MSSM-weak[%] 


-1.22 


-1.57 


-1.60 



Table 2: LO cross section in the improved Born approximation ctiba as defined in (j2.53p . 
as well as corrections 5 relative to cjiba, for pp —>■ (bb) h'^/H'^/A°-|-X at the LHC {^/s = 
14 TeV) in the SPS 4 scenario. The MRST2004qed PDFs and NNLO-QCD running 
for TTib have been adopted, the renormalization and factorization scales have been set to 
fiR = M^o and fip = M^o/A, respectively, and tp has been renormalized in the DR scheme. 
"QCD" denotes the NLO QCD corrections only, "QED" denotes all photonic corrections 
only, and "MSSM-QCD" and "MSSM-weak" comprise only the QCD and weak effects in 
the MSSM, respectively, that remain after absorbing the large-t^ effects in the LO cross 
section. 



the Higgs masses, in contrast to the default setting in FeynHiggs. Including the imaginary parts 
has no visible effect on M^o and shifts Mjjo by only approximately 200 MeV.) 

3.4 The MSSM cross sections 

Within the MSSM, let us first focus on the radiative corrections and total cross sections in the 
SPS 4 benchmark scenario using the DR scheme for the renormalization of tp. 

In Table [2] we present the relative radiative corrections 5 defined with respect to the im- 
proved Born approximation o"iba defined in Eq. (I2.53p . The full cross-section prediction includ- 
ing summations and the remaining non- universal 0{as) and 0{a) corrections is thus given by 
cr = aiBA X (1 (5qcd + Sqeb +<^mssm-qcd + ^MSSM-weak)- Note that we have removed cor- 
rections from our full calculation that are taken into account through the use of Ocfr in ctiba to 
avoid double counting. As can be seen from Table [21 the bulk of the MSSM-QCD and -weak cor- 
rections can indeed be absorbed into the above definition of ctiba- The remaining non- universal 
corrections (5mssm = <5mssm-QCD + '5MSSM-weak in the complete MSSM calculation turn out to be 
quite small, below approximately 2%. We note, however, that the size of the corrections (^mssm 
depends quite sensitively on the numerical value of the input b-mass. We shall discuss this point 
in more detail below. 

In Table [3l we show the cumulating effect of the various higher-order corrections on the effec- 
tive b-mass (as defined after summation of t/3-enhanced terms, mh^mi,{l — A\j/{tpta))/{l + Ab) 
for h*^ production, etc.) and the resulting cross sections. As already observed within the SM, 
the QED corrections are generally very small after mass factorization. The summation of the 
t^-enhanced MSSM-QCD and MSSM-weak corrections, encoded in and A^^^^, respectively, 
has a significant effect on the cross sections for H'' and A*^ production. The light Higgs boson h", 
on the other hand, is SM-like in the SPS 4 scenario and the summation of terms oc tp has thus 
no sizeable impact on the cross section. Employing a loop-improved effective mixing angle agg is 
numerically relevant only for h*^ production b GcausG Sq^ ^ — Xjtj^ is small and. Gvcn a small shift 
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hO 


HO 


AO 


mb[GeV] 


a[ph] 


mb[GeV] 


a[ph] 


mb[GeV] 


a[ph] 


QCD 


2.80 


0.97 


2.55 


24.12 


2.55 


24.13 


+QED 


2.80 


0.97 


2.55 


24.07 


2.55 


24.09 




2.72 


0.92 


1.95 


14.14 


1.95 


14.15 


_^^weak 


2.75 


0.94 


2.24 


18.66 


2.24 


18.67 


+ sin(ae//) 


2.75 


0.88 


2.24 


18.66 


2.24 


18.67 


full calculation 


2.75 


0.87 


2.24 


18.43 


2.24 


18.44 



Table 3: The NLO MSSM cross section pp ^ (bb) IiVhVA^ + X at the LHC {^/s = 
14 TeV) in the SPS 4 scenario including the cumulative corrections due to the different 
classes of corrections (PDFs, scale setting, and renormalization as in Table [2]). We 
also quote the effective bottom mass entering the respective effective Yukawa coupling, to 
quantify the impact of the summation. "QCD" denotes the NLO QCD corrections only, 
"QED" denotes the addition of all photonic corrections, and A^'''^'^ refer to the effect 
of summing enhanced terms in SUSY-QCD and the weak MSSM, respectively, sin(Q;eff) 
denotes the additional improvement of the Born cross section by using the effective mixing 
angle Ocs- Finally, we give the result for the full calculation in the MSSM where all 0{as) 
and 0{a) corrections are included. 



a Oeff has a sizeable effect on Sq;. The cross sections in the last-but-one row of Table [3] cor- 
respond to the improved Born approximation dressed with QCD and QED corrections. The full 
MSSM cross section including all summations and the remaining non-universal 0{as) and 0{a) 
corrections is displayed in the last row of the table. As discussed previously, the non-universal 
supersymmetric corrections turn out to be small at a level of a few percent. Note that the b- 
mass values in the last row of Table [3] have been used to calculate the cross sections tabulated 
in Table O (While mb = 2.75 GeV and mb = 2.24 GeV have been used to calculate ctiba for 
h*^ and if/A^ production, respectively, mb = 2.24 GeV has been used as input for all relative 
corrections, as discussed before.) 

In Fig. [2] we show the impact of the complete supersymmetric 0{as) and 0{a) corrections, 
^MSSM = (^MSSM-QCD + ^MSSM-weak, defined relative to the improved Born approximation cjiba 
(j2.53p . for different values of the on-shell mass M^o. All other MSSM parameters are kept 
fixed to the SPS 4 values. The size of the non-universal corrections does not exceed 3% for 
H''/A'^ production except for special model parameters, where the Higgs masses are close to 
the production threshold for pairs of sparticles. The peaks in the corrections correspond to 
neutralino, chargino, or sfermion thresholds. Thresholds which are too narrow, however, are not 
displayed in Fig. [2j These unphysical singularities can be removed by taking into account the 
finite widths of the unstable sparticles (see e.g. Ref. [39]). Note that the peaks for h*^ production 
are induced by the finite parts of the counterterms in (|2.34p . This proliferation of unphysical 
singularities can be avoided by choosing the DR scheme for the A'' wave function renormalization, 
as is default in FeynHiggs. 

For small M\o when the masses of the neutral Higgs bosons are almost degenerate, the 
effects from the loop-induced Higgs mixing become extremely large. To go beyond the effective 
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200 300 400 500 600 700 800 900 1000 

MAo[GeV] 

Figure 2: Full MSSM corrections (^mssm = <Jmssm-qcd + f^MSSM-woak defined relative 
to (TiBA as a function of the Mpfi pole mass in the DR scheme for tfj. All other MSSM 
parameters are fixed to their SPS 4 values. 
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<5mssm [%] 





1 hO 
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200 300 400 500 600 700 800 900 1000 

MAo[GeV] 

Figure 3: Full MSSM corrections (5mssm = (^mssm-qcd + (^MSSM-weak as in Fig. [2l 
However, here a^s (!2.52p in ctiba is calculated from self-energies at = 0. 
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mixing angle approximation in this region, one would have to include corrections to the ofF-shell 
bbh'^/bbH" vertex for H^'/h*^ production, respectively. We will not address this issue, which is 
part of a two-loop calculation, in this work. Therefore, we truncate Fig. [2] at M^o = 150 GeV. 

For very large M^o the relative corrections to h*^ production increase to up to 10%. However, 
in this parameter region, the size of 5mssm for h*^ depends very sensitively on the definition of 
the effective mixing angle a^s employed in (Tiba- As defined in (j2.52p . a^s inherits a distinct 
dependence on M^o because some of the self-energies are evaluated at momentum scales of the 
order of Mjjo ~ M^o. This dependence is not present in the complete result and, thus, has 
to be compensated by (5mssm- Evaluating the self-energies that enter a^s (j2.52p at = 0, 
the peak structure for large M^o in Fig. [2] is absent, and the size of the correction is 2-3% for 
300 GeV < M^o ^ 1000 GeV, as shown in Fig. [3l However, this approximation breaks down 
at small values of M^o where the corrections due to Higgs mixing become large. Note that in 
any case h^ is SM-like at large M^o so that h^ production is most likely of no phenomenological 
relevance. 

The two-loop improvement of the self-energies of the CP-even Higgs bosons, as contained in 
the Z factors of (j2.4ip and (j2.42p . has a negligible effect on Jmssm for H*^ and A'' production. 
Only for h*^ production and M^o ^ 350 GeV the difference between a one-loop and two-loop 
treatment is larger than 1% and can reach the 5% level for M^o ~ 200 GeV. However, the 
two-loop improvement plays an important role for the precise determination of the Higgs-boson 
pole masses entering the kinematics of the process [40] . 

It is important to emphasize that the size of the non-universal MSSM one-loop corrections 
^MSSM at large tjs is quite sensitive to the choice of the bottom-mass input value. This is caused 
by terms that grow as t^p but are not included in the summation of t^-enhanced terms. For the 
SPS 4 scenario, the b-mass input sensitivity is shown in Fig. [H In the DR scheme for (upper 
panel of Fig. S]) the absolute size of the non-universal corrections varies between approximately 
zero and —6% for YL^/A^ production, depending on whether a massless approximation, the 
running mass or the pole mass is chosen as b-mass input. The sensitivity is even larger for the 
light Higgs h*^. However, for h*^ production (5mssm depends very sensitively on the definition of 
OefT) as mentioned before. Although we assume that the running mass including the summation 
of t/3-enhanced terms is a sensible choice, the sensitivity of the NLO correction to the b-mass 
input constitutes a theoretical uncertainty which can not be resolved at the NLO level. 

Comparing the DR and the DCPR schemes for the renormalization of tp, the size of the 
corrections is accidentally very similar if the running b-mass is chosen as an input. The same 
observation holds for the Higgs masses, which differ by less than 500 MeV. However, because 
of the large mb dependence the corrections in the different renormalization schemes can differ 
significantly in general. This can be seen by comparing the upper and lower panels of Fig. [H For 
the Higgs masses we also find larger differences up to 4 GeV for m-b = 4.2 GeV. Of course, the large 
scheme difference of the relative correction is partially compensated by the corresponding change 
in the DCPR input value for to be calculated from the DR value tp = 50. In the massless-b 
approximation, we find = 51.78, while = 47.00 for mb = 4.2 GeV. The resulting scheme 
dependence of the total cross-section prediction for H*^ and A*^ production is thus moderate 
(below 1%) even for large mb. For h*^ and large mb, the residual scheme dependence can reach 
up to 3%. To compare the cross sections, we have used the on-shell mass for A*^ as computed in 
the DR scheme as input in both schemes. 

We have verified that the size of (5mssm for H*^ or A*^ production does not change significantly 
when the input values of the soft breaking parameters for the sfermions or the soft breaking 
parameters for the gauginos are varied around their SPS 4 values by up to a factor of 2, unless 
some sparticles become unreasonably light. 
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Figure 4: Full MSSM corrections (5mssm = <Jmssm-qcd + (^MSSM-wcak defined relative to 
(TiBA as a function of the rrib input in the DR scheme (upper panel) and DCPR scheme 
(lower panel) for i^. The corrections for H'^ and A'' lie on top of each other. 
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Results for the alternative benchmark scenario SPS lb with tp = 30 are presented in Ap- 
pendix [Bl The conclusions reached for the SPS 4 scenario essentially hold for SPS lb as well. 

4 Conclusions 

We have performed a complete calculation of the 0{as) and 0{a) corrections to associated 
bottom-Higgs production through bb fusion in the MSSM. This next-to-leading order prediction 
is improved by including two-loop corrections, as provided by FeynHiggs, and by the known 
summation of f^-enhanced corrections. We have presented details of the calculation and dis- 
cussed numerical results for MSSM Higgs-boson production at the LHC in two supersymmetric 
scenarios. The leading supersymmetric higher-order corrections, in particular the t^-enhanced 
contributions, can be taken into account by an appropriate definition of the couplings and the 
running b-mass in an improved Born approximation. The remaining non- universal corrections 
are small, typically of the order of a few percent. The quality of an improved Born approxima- 
tion, however, can only be judged with a full 0{as) and 0{a) calculation. Although we assume 
that the running mass including the summation of t^-enhanced terms is a sensible choice for the 
input b-mass, alternative choices can lead to non-universal corrections as large as 10%. 

With the results presented, the impact of the complete MSSM corrections to neutral Higgs 
boson decays into bottom quarks might also be updated. 

Our results show that the difference between a properly defined improved Born approximation 
and the complete NLO calculation, which is improved by leading NNLO effects, is smaller than 
other theoretical uncertainties resulting from residual scale dependences, errors on the b-quark 
mass, and parton distribution functions. 
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A SPS benchmark scenarios 

For the SPS benchmark [38] scenarios discussed in this work we use the following input for t^, 
the mass of the CP-odd Higgs boson Mao, the supersymmetric Higgs mass parameter /i, the 
electroweak gaugino mass parameters Mi^2, the gluino mass rrig, the trilinear couplings ^T,t,b) 
the scale, at which the DR input values are defined /i/j(DR), the soft SUSY-breaking parameters 
in the diagonal entries of the squark and slepton mass matrices of the first and second generation 
Mfi (where i = L,R refers to the left- and right-handed sfermions, / = q,l to quarks and 
leptons, and f = u,d,e to up and down quarks and electrons, respectively), and the analogous 
soft SUSY-breaking parameters for the third generation M?p: 
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SPS 4 


SPS lb 






SPS 4 


SPS lb 




50.0 


30.0 




MgLlGeY] 


732.2 


836.2 


MAo[GeV] 


404.4 


525.5 






713.9 


803.9 


fi[Ge\] 


377.0 


495.6 




MuRlGeY] 


716.0 


807.5 


Mi[GeV] 


120.8 


162.8 




MiLiGeY] 


445.9 


334.0 


M2[GeV] 


233.2 


310.9 




Mei?[GeV] 


414.2 


248.3 


mg[GeV] 


721.0 


916.1 






640.1 


762.5 


^^[GeV] 


-102.3 


-195.8 




M|f[GeV] 


673.4 


780.3 


At[Ge\] 


-552.2 


-729.3 






556.8 


670.7 


Ab[GeV] 


-729.5 


-987.4 




M3f[GeV] 


394.7 


323.8 


/xij(DR)[GeV] 


571.3 


706.9 






289.5 


218.6 



B Results for SPS lb 

In the SPS lb scenario the two-loop Higgs masses are given by M^o = 117.67 GeV, Mjjo = 
525.69 GeV, and M^o = 525.66 GeV in the DR scheme for t^. Here, the masses in the DCPR 
scheme also differ by not more than 500 MeV and the discrepancy does not increase as drastically 
as in SPS 4 when the input value for the b-mass is changed with respect to its default value 
rrib = 2.30 GeV. As can be seen from Table [H the non-universal corrections in the MSSM for 
A'^/H'^-production are even smaller than for SPS 4. Because the masses of A^/H*^ are larger and 
tj3 = 30 is smaller than in SPS 4, the total production cross sections are also smaller. The cross 
sections including the cumulating effect of the different higher-order corrections are shown in 
Table [5l The generic structure of (5mssm as a function of Af^o with all other SPS lb parameters 
fixed. Fig. El does not differ from SPS 4. 







HO 


A" 


0"IBA[pb] 


0.59 


1.88 


1.88 


^qcd[%] 


35.98 


19.02 


19.02 


Sqed[%] 


-0.13 


-0.26 


-0.26 


(^mssm-qcd[%] 


-0.06 


0.11 


0.11 


'^MSSM-weak[%] 


2.73 


-0.35 


-0.35 



Table 4: LO cross section in the improved Born approximation (Jib a as defined in (j2.53p . 
as well as corrections 6 relative to ctibA) for pp (hh)]i^ /H^ / +X at the LHC (t/s = 
14 TeV) in the SPS lb scenario. See Table [2] for details. 
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mb[GeV] 


a[pb] 


mb[GeV] 


a[pb] 


QCD 


2.79 


0.85 


2.50 


2.64 


2.50 


2.64 


+QED 


2.79 


0.84 


2.50 


2.63 


2.50 


2.63 




2.76 


0.83 


2.08 


1.82 


2.08 


1.82 


1 A weak 


2.77 


0.83 


2.30 


2.24 


2.30 


2.24 


+ sin(ae//) 


2.77 


0.80 


2.30 


2.24 


2.30 


2.24 


full calculation 


2.77 


0.81 


2.30 


2.23 


2.30 


2.23 



Table 5: The NLO MSSM cross section pp ^ (bb) hO/RO/AO + X at the LHC {^/s = 
14 TeV) in the SPS lb scenario. See Table [3] for details. 
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Figure 5: Full MSSM corrections (5mssm = <Jmssm-qcd + (^MSSM-weak defined relative 
to (TiBA as a function of the M^o pole mass in the DR scheme for tj3. All other MSSM 
parameters are fixed to their SPS lb values. 
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